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ABSTRACT 

The exchange of meteorites among the terrestrial planets of our Solar System is a well 
established phenomenon that has triggered discussion of lithopanspermia within the Solar 
System. Similarly, could solid material be transferred across planetary systems? To address 
this question, we explore the dynamics of the transfer of small bodies between planetary 
systems. In particular, we examine a dynamical process that yields very low escape veloci- 
ties using nearly parabolic trajectories, and the reverse process that allows for low velocity 
capture. These processes are chaotic and provide a mechanism for minimal energy transfer 
that yield an increased transfer probability compared to that of previously studied mecha- 
nisms that have invoked hyperbolic trajectories. We estimate the transfer probability in a 
stellar cluster as a function of stellar mass and cluster size. We find that significant amounts 
of solid material could potentially have been transferred from the early Solar System to our 
nearest neighbor stars. While this low velocity mechanism improves the odds for interstellar 
lithopanspermia, the exchange of biologically active materials across stellar systems depends 
greatly upon the highly uncertain viability of organisms over the timescales for transfer, 
typically millions of years. 
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1. INTRODUCTION 



From the collection of thousands of meteorites found on Earth, there are about 20 
that have been identified as having a Martian origin, and a similar number that originated 
from the Moon. The study of the dynamical evolution of these meteorites agrees well with 
the cosmic ray exposure time and with the frequency of landing. No meteorites were ever 
found on the Moon by the Apollo mission. In 2005, the rover Opportunity encountered in 
Mars the first meteorite on another Solar System body, identified as an iron-nickel mete- 
orite. These findings, together with dynamical simulations (Gladman 1997; Dones et al. 
1999; Mileikowsky et al. 2000), indicate that meteorites are exchanged among the terrestrial 
planets of our Solar System at a measurable level. Because sufficiently large rocks may pro- 
tect dormant microorganisms from cosmic ray exposure and from the hazards of the impact 
at landing, it has been suggested that the exchange of microorganisms living inside rocks 
could take place among the Solar System planets, a phenomenon known as lithopanspermia. 
In fact, new laboratory experiments have confirmed that several microorganisms (bacterial 
spores, cyanobacteria and lichen) embedded in martian-like rocks could survive under shock 
pressures similar to those suffered by martian meteorites upon impact ejection (Stoffler et 
al. 2007; Horneck et al. 2008). Under this scenario, life on Earth could potentially spread to 
other moons and planets within our Solar System, and/or life on Earth could have an origin 
elsewhere in our Solar System. 

Melosh (2003) investigated the probability of lithopanspermia taking place amongst the 
stars in the solar local neighborhood. He found that even though numerical simulations show 
that up to one-third of all the meteorites originating from the terrestrial planets are ejected 
out of the Solar System by gravitational encounters with Jupiter and Saturn, the probability 
of landing on a terrestrial planet of a neighboring planetary system is extremely low because 
of the high relative velocities of the stars and the low stellar densities. He concluded that 
lithopanspermia among the current solar neighbors is "overwhelmingly unlikely" . 

In a subsequent paper, Adams & Spergel (2005) pointed out that the majority of stars, 
including the Suiil], are born in stellar clusters with N = 100-1000 members and a typical 
radius of R = lpc{N/100y/'^ (Lada & Lada 2003 and Carpenter 2000). In such an envi- 
ronment, the probability of transfer would be higher due the larger stellar densities and 
smaller stellar relative velocities compared to those for field stars (and for the current solar 



^In the Solar System, the existence of short-lived isotopes (^°Fe and ^^Al) in primitive meteorites has 
been interpreted as indication of nearby a supernova explosion shortly before the Solar System solids started 
to accrete. This has been considered as evidence that the Sun was born in a cluster environment. 
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neighborhood) . The timescale for planet formation and the dispersal time of the clusters are 
comparable (10-100 Myr); therefore, it could be possible that solid material be transferred 
before the cluster disperses. Adam & Spergel (2005) estimated the probability of transfer 
of biologically-active remnants between planetary systems within a cluster by using Monte 
Carlo simulations, assuming that the stars are mostly in binary systems (which increases the 
cross-section). Adopting typical ejection speeds of ~5 km/s, they found that the expected 
number of successful lithopanspermia events per cluster is ~ 10~^; for lower ejection speeds, 
~ 2 km/s, this number is 1-2. 

Because there is a significant increase in the number of possible lithopanspermia events 
with decreased ejection velocity, it is of interest to study a very low energy mechanism 
with velocities significantly smaller than those considered in Adams & Spergel (2005). This 
mechanism was described by Belbruno (2004) in the mathematical context of a class of 
nearly parabohc trajectories in the restricted three-body problem. The escape velocities of 
these parabolic-type trajectories are very low, ~ 0.1 km/s, substantially smaller than the 
mean relative velocity of stars in the cluster, and the remnant escapes the planetary system 
by slowly meandering away. This process of "weak escape" is chaotic in nature and its 
study requires the use of methods of chaos theory, lying beyond the reach of Monte Carlo 
simulations employed in previous studies. 

"Weak escape" is a transitional motion between capture and escape. For it to occur, 
the trajectory of the remnant must pass near the largest planet in the system. "Weak 
capture" is the reverse process, when a remnant can get captured with low velocity by 
another planetary system. The fact that the escape velocities of the remnants we consider 
here are small, enhances the probability that a remnant can be weakly captured by another 
planetary system due to lower approach velocities of the stars encountered. 

In § 2, § 3 and § 4, we describe the mathematics of weak escape and weak capture (based 
on Belbruno 2004), and in § 5 we disscuss its astronomical application to the study of the 
slow chaotic transfer of solid material (remnants) between planetary systems within a star 
cluster: § 5.1 describes the location of the region supporting weak capture and weak escape 
for a star of a given mass; §5.2 constraints the range of stellar masses that allow minimal 
energy transfer (weak transfer) to take place; §5.3 estimates the probability of weak capture 
by a star of a given mass; and §5.4 calculates the number of weak transfer events. The focus 
of the paper is the study of the transfer of solid material between planetary systems, without 
regard to whether or not they might contain material of biological interest. Our motivation, 
however, is to shed light on the following questions: could the building blocks of hfe on 
Earth have been transferred to other planetary systems within the first 10-100 Myr of the 
Solar System evolution, when the Sun was still embedded in its maternal aggregate?; and 
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vice versa, could life on the Solar System have been originated beyond its boundaries? This 
is addressed in §5.5, where we apply the results in the previous sections to lithopanspermia, 
focusing in particular on solar-type stars. 



2. MODEL 



We first describe the mathematics of weak escape and weak captm^c. Readers interested 
in the astronomical application only could proceed to §5. We define a general planetary 
system, S, consisting of a central star. Pi, and a system of N planets, Pj,i = 2,...,N 
(N > 3) on co-planar orbits that are approximately circular. Their labeling is not reflective 
of their relative distances from Pi. We assume that the mass of the star, mi, is much 
larger than the masses of any of the planets, rrii (mi 3> m,,? = 2,...,N), and that the 
mass of one of the planets, P2, is much larger than the sum of the masses of all the other 
planets, m2 3> mi,i = 3, N (this condition is fulfilled in the case of our Solar System 
with P2 = Jupiter). These assumptions reduce S* to a much simpler system without the 
loss of generality. We consider a remnant, Pq, whose mass, mo, is negligible with respect 
to Pj, i — 1,...,N; we assume that Pq orbits the star in the same plane as the planets 
without affecting their orbits, and that it spends most of the time far from Pi. Because 
m2 ^ "rriiji — 3,...,N, the gravitational perturbation of P2 on the motion of Pq is the 
dominant one, and we will ignore the perturbations due to the other planets. This reduces 
the motion of Pq to that of a three-body problem between Po,Pi and P2, where Pi and P2 
are moving in approximately circular orbits about their common center of mass. Because 
mi ^ m2, we can view Pi as fixed, with P2 orbiting around it in a circular orbit at constant 
radial distance A {lAU < A < 500 AU) and orbital frequency uj. This defines the classical 
"planar circular restricted three-body problem" for the motion of Pq. 

The differential equations for the restricted three-body problem are well known. We 
write them in dimensionless form; the details can be found in Belbruno (2004). Without 
loss of generality, we choose a reference frame with origin at the center of mass of the 
Pi,P2 system; we choose units such that A = 1, u = 1, mi = 1 — /i and m2 = ^, where 
fi = 1712/ {mi + m2) > and /i ^ 1. Under these assumptions, the period of motion of P2 
around Pi is T = 27r. We will refer to this as the reduced Solar System, Sq. In our reduced 
Solar System, Pi is the Sun, P2 is Jupiter and ji ^ 0.001. In inertial coordinates {Qi,Q2) 
the differential equations for the motion of Pq can be written as 



Q — ^Q, 



(1) 
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where we have used the notation ' = ^, and Q = {Qi,Q2), assumed to be a vector, and 
= fg, with 

n = i^ + A (2) 

ri r2 

ri(t) = V(Qi-H2 + (Q2-/is)2, (3) 

r2{t) = V(Qi + (l-/i)c)2 + (Q2 + (l-/x)s)2, (4) 

where ri, r2 are the distances of Pq to Pi and P2, respectively, c = cos(t), s = sin(t), and the 
position of Pi is given by /i(c, s) and that of P2 is given by —(1 — /i)(c, s). 

These differential equations can also be writen in a barycentric rotating coordinate 
system {xi,X2) with orbital frequency u. In this case, they form an autonomous system with 
an energy integral, the Jacobi energy, J = J{x,x), where x = {xi,X2) is a vector and J is 
a function on the four-dimensional phase space {x,x). Along a solution, (x(t),x(t)), J is a 
constant of the motion and defines a three-dimensional surface constraining the motion of 
Pq called the Jacobi surface. In the rotating system. Pi and P2 are fixed and it is convenient 
to place them on the Xi-axis, with Pi at xi = fi and P2 at xi = — 1 + /i, in which case, 

J=_|i;|2+|x|2 + /x(l-/x)+2fi. (5) 

where |x| represents the standard Euclidian norm of x, and the additive term — yu) is 
present so that the values of the Jacobi energy, J = C, are normalized. 



3. WEAK CAPTURE AND ESCAPE 



In this section, we introduce the concepts of weak capture and escape (for a detailed 
discussion we refer to Belbruno (2004, 2007b). A convenient way to define the capture of 
Pq with respect to Pi or P2 is by using the concept of "ballistic capture". We define the 
two-body Kepler energy, E^, of Pq with respect to one of the bodies Pk, k = 1,2: 

Ek = -v 6 

2 Tk 

where v is the velocity of Pq relative to Pk- Ballistic capture takes place when < 0. More 
precisely, let </)(t) = {Qit),Q{t)) be a solution of Eq.([l]) for Pq, and assume no collisions take 
place, i.e. > 0. Pq is ballistically captured by Pk at time t = t* if Ek{4>(t*)) < 0; similarly, 
Pq escapes P^ at time t = t2, if there exists a finite time interval [^1,^2], with t2 > ti where 
Ek{(f){ti)) < and Ek{(f){t2)) > 0. For notation, we set Ek{(f){t)) = Ek{t). 
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Now we consider the motion of Pq around P2 (this is known as Hill's problem). If the 
motion of Pq is initially circular and occurs sufficiently close to P2, it is found that the motion 
is generally stable and E2 remains negative for all times (this was proven by Kummer(1983) 
applying the Kolmogorov-Arnold-Moscr theorem). However, if Pq is not sufficiently near 
to P2, then the stability will break down and the motion will substantially deviate from a 
circular orbit, leading to escape from P2 due to the gravitational perturbation of Pi. 

Of particular interest is temporary ballistic capture. This occurs when Pq is ballistically 
captured by P2 over the finite interval [ti,t2]j i-e., when E2{t) < 0, for ti < t < t2, ^2(^1) = 
0,i?2(^2) = 0, and E2{t) > for t < ti and t > t2. We say that Po is "pscudo-ballistically 
captured" at a time t = t* when E2{t*) ^ 0, i.e. where the Kepler energy is slightly 
hyperbolic. 

Because temporary ballistic capture with respect to P2 is generally unstable and chaotic 
in nature (Belbruno 2004; Garcia & Gomez 2007) and lies at the transition between capture 
and escape, it is referred to as "weak capture" . For weak capture to occur, the Jacobi energy 
C must be small enough so that the velocity of Pq is sufficiently large at a given distance 
from P2 to be in this transition state. 

As is evident by numerical integration of Pq about P2, weak capture generally occurs 
for relatively short time spans. At = t2 — ti. For example, in the case of the Earth- Moon 
system (// = 0.012), weak capture about the Moon occurs for time spans of days or weeks; 
in the case of the Sun- Jupiter system, the timescale is months or a few years (Belbruno & 
Marsden 1997; Belbruno 2007a) 

A numerical approach to estimate the region around the Moon supporting weak capture 
was developed by Belbruno (2004). This region was termed the weak stability boundary, 
W, and it is the location in the phase space where the motion of Pq with respect to P2 
lies between capture and escape. Recent results by Garcia & Gomez (2007) have precisely 
determined W using the numerical approach and have shown that it is closely related to 
a complicated fractal region in phase space supporting unstable chaotic motion. However, 
the location of W can be can be approximately estimated analytically in a straight forward 
manner by a set W described as follows: 

Consider the barycentric rotating coordinates x = {xi,X2). To be in weak capture, Pq must 
have a suitably large velocity magnitude, \x\, that depends on the velocity direction at a 
given distance r2 from P2. This is equivalent to the requirement that the Jacobi constant 
C needs to be in a suitably small range: we require that C ~ Ci, where Cj is the value of 
J at the classical Lagrange points Lj (which are the locations at which = 0). We label 
the three coUinear Lagrange points from left to right (in the rotating barycentric reference 
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frame defined in § 2) with index 1,2, and 3; then the values of the Jacobi energy, Cj , at each 
of these points has the ordering: C2 > Ci = 3. As C decreases, the allowed range of motion, 
X, of Pq expands. When C goes just below Ci, the Hill's regior|^ around P2 opens near the 
position of Li so that Pq can move out of the Hill's region around P2, and infinitely far away 
from both Pi and P2 (see Figure [1]- i.e. Ci is the minimal Jacobi energy for escape). If 
we decrease C further, C < 3 < Ci, then the zero velocity curves disappear and the Hill's 
region becomes the entire a;iX2-plane. 





Fig. 1. — Basic Hill's regions (in white): Starting from left to right and top to bottom, C 
has the values: C > C2; Ci < C ~ C2; C ~ d; C3 < C < Ci and C < 3. 

As is described in Belbruno (2004), the set W has two components. We and Wh, as 
follows. 

W = WeUWh, (7) 
We = {(x,x) e M^|E2 < 0, J = < C < Ci,f2 = 0}, (8) 
Wh = {{x,x) eR^\E2^0,J = C,C* <C <Ci}, (9) 

where the constant C* is determined so that W exists and depends on the value of fi. Weak 
capture occurs on We (associated with osculating elliptic orbits about P2), and pseudo- 
ballistic capture (equivalently pseudo-ballistic escape) occurs on Wh (associated with slightly 
hyperbolic orbits). 

For weak capture we need to consider the set We, 

We = J{C) n S n a, (10) 



^The Hill's region for a given value of C is the projection of the Jacobi surface onto the Xi^xi space, 
which yields locations where Pq is allowed to move. 
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where J = {{x,x) e R V = C}, E = {{x,x) E ^."^{£2 < 0}, a = {(x,x) e R^lh = 0}. We 
is equivalent to a set of osculating elliptic orbits about P2 with osculating pericenter q2, in 
a two-dimensional annular region described by 

g2 = /(^2,e2), (11) 

where 62 {0 ^ O2 < 27t) is the polar angle measured with respect to the Pi,P2 axis) in a 
P2-ccntcred rotating coordinate system, the function / is periodic in 62 of period 27r, and 62 
is the osculating eccentricity of the orbit of Pq with respect to P2 (0 < 62 < 1). For example, 
in the case where C ~ Ci and q2 ~ 0,n ~ 0, this relationship takes the form. 



[1- 62)1^3 
3t 



.2 - 'ir^. (12) 



3' 



The set Wh plays a key role on the existence of chaos, as it is proven that a hyperbolic 
invariant set associated to parabolic motion with respect to Pi exists on the set Wh, giving 
rise to chaotic dynamics (Belbruno 2004, pages 186-191). 

A more accurate estimation of W has recently been obtained for a large range of the Ja- 
cobi energies by Belbruno et al. (2007b) through the visualization of special two-dimensional 
Poincarc sections, generally showing resonant tori lying within a chaotic sea. These sections 
yield a more accurate representation of W, suggesting that it may consist of transverse 
homoclinic points, related to the intersections of the invariant manifolds of the Lyapunov 
periodic orbits associated to the unstable Lagrange points. It is demonstrated that during 
weak capture, Pq moves around Pi while transitioning between various resonant states with 
respect to P2. However, these results apply to elliptic- like resonant motions that are not 
sufficiently energetic for the purposes of this paper. In the following, we will be considering 
the parabohc motions associated with Wh- 

In our Solar System, the existence of the weak stability boundary and the viability 
of weak capture was demonstrated in practice by the Japanese spacecraft Hiten: using 
a trajectory designed by Belbruno (1990, 1993, 2007a), Hiten was captured into an orbit 
around the Moon in 1991 without the use of rockets to slow down. Weak capture at the 
Moon was also achieved in 2004 by the ESA spacecraft SMARTl (Racca 2003, Belbruno 
2007a). In both cases, the capture was chaotic and unstable because the orbit lay in the 
transition between capture and escape. In another application, weak escape from the Earth's 
L4 (or L5) was invoked to suggest a low energy transfer to the Earth for the hypothetical 
Mars-sized impactor that is thought to have triggered the "giant impact" origin of the Moon 
(Belbruno & Gott 2005). 
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4. PARABOLIC MOTION AND CHAOS 

In this section, we describe the mechanism for low velocity escape from 5* that is chaotic 
and applies to the set of parabolic trajectories around Pi. The reverse process yields low 
velocity chaotic capture into S. Consider a solution Q(t) of ([T]) for Pq in the case of the 
reduced Solar System So for the restricted three-body problem: Pq is on a parabolic trajectory 
with respect to Pi when limt^±oo \Q\ = oo and limf^-too |<5| =0. In the absence of any of the 
planets, the motion of Pq is that of a standard parabola around Pi, and the Jacobi energy 
has the value C = ±2 a/2, where + is for retrograde parabolic trajectories, and - is for direct 
trajectories. 

In the case of the Solar System, S, because the masses of the planets are very small with 
respect to that of Pi, the motion of Pq will be slightly perturbed and the trajectories that lie 
very close to parabolic trajectories will be able to escape Pi with very small velocities; we refer 
to these as "pseudo-parabolic trajectories". These trajectories satisfy \Q\ R, \Q\ a, as 
t — > Tr, and, without loss of generality, we assume that Pq started at periapsis with respect 
to Pi at time t = and traveled a distance P in a sufficiently long time Tr > 0. The value 
of a is small and coincides with the escape velocity with respect to Pi at the distance R, 
a = ^j2Gmil R\ R is chosen sufficiently large so that Pq can escape S when it has a velocity 
close to cr. 

For notational purposes, the term "parabolic trajectory" refers to the solution of ([1]) 
for 5*0 in scaled dimensionless coordinates and the precise definition of parabolic motion; 
whereas, when using the term "pseudo-parabolic trajectory", we are considering system S 
in unsealed dimensional coordinates. We now describe the main results for 5*0, followed by 
the results for S. 

Main Results for So'- 

Let (j)(t) = {Qit),Q{t)) be a parabolic solution of ([T]) for Pq in a system where = 
and where Pq does not collide with Pi (i.e. ri > 0). Consider the set of parabolic trajectories 
with \C\ ~ 2^2 ^ 2.83, that pass through the location of P2 (now with zero mass because 
/i = 0). Belbruno (2004) found that for < ;U <^ 1, if Pq passes sufficiently close to P2 
without collision (r2 > 0) and slightly beyond the distance A, then: (a) when Pq passes 
close to P2, it is slightly hyperbolic and in the set Wh (i.e. pseudo-weakly captured) around 
P2; (b) Pq never collides with Pi or P2; (c) the periapsis location of Pq with respect to P2 
is approximately the periapsis location with respect to Pi; (d) there exists a set of positive 
measure of parabolic trajectories (i.e. parabolic trajectories exist); and (e) the motion of Pq 
is chaotic. 
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There are two possible outcomes: 

1. Pq parabolically escapes Pi, i.e. — > oo and jQl — > as i — > oo. 

2. Pq does not escape but instead, after traveling a large finite distance, falls back towards 
Pi, passing again by P2 in a slightly hyperbolic orbit w.r.t P2 in which Pq is pseudo- 
weakly captured by P2, after which it flies around Pi and out again toward \Q\ — 00. 

In the second case, every time Pq falls back to fly by P2, the argument of periapsis of 
the pseudo-parabola with respect to Pi is chaotic in nature and can take on a random value. 
If Pq were to keep falling back to P2 again and again, it would be permanently captured 
by Pi. However, it can be proven that this docs not happen (in mathematical terms this is 
because the set of orbits leading to permanent capture are of zero measure); on the contrary, 
escape eventually takes place (i.e. there is a set of positive measure of parabolic trajectories 
that will escape to infinity after flying by P2). In addition, in this case, there is a constraint 
on the fly- by distance q2 as a function of If we consider the case of a direct fly by, the 
velocity of Pq with respect to Pi when it flies by P2 is V ^ 1.414 (in dimensionless units) - 
this is because the fly-by is close to P2 and therefore has a distance A 1 to Pi, and hence, 
V ~ y^2/A ?a 1.414); therefore, the velocity of Pq with respect to P2 is approximately 0.414, 
and the Kepler energy of Pq with respect to P2 is E2 ~ 0.08 — {/j/r2). A slight hyperbolic 
fly-by implies that E2 ~ 0, or equivalently, r2 ~ In the case of the Solar System, if P2 
is Jupiter (/x = 0.001) then r2 - 0.013 = 10,172,800 km, while for Neptune (/x = 0.00005), 
rs ~ 0.0006 = 2,692,800 km. 

Main Result for S: 

The parabolic motion previously described for Pq in the system Sq can be translated to 
the pseudo-parabolic motion of Pq in the system S in the following manner. There is a set of 
positive measure of pseudo-parabolic trajectories with respect to Pi that pass by P2 without 
colliding with it, and where Pq is slightly hyperbohc with respect to P2. The motion of such 
trajectories is sensitive due to the fly- by. As Pq approaches P2, it achieves a periapsis with 
respect to P2 slightly beyond P2's orbit, with energy at periapsis slightly hyperbolic with 
respect to P2. At that location, it is also approximately at the periapsis with respect to 
Pi. Pq then flies outward away from Pi for a long period of time Tr, reaching a distance R 
with a velocity ~ a. Pq pseudo-parabolically escapes Pi with approximate parabolic escape 
velocity magnitude a, when ri increases beyond R as t increases beyond Tr. 

Weak Stability Boundary Around Pi 
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The pseudo-parabolic escape of Pq occurs for R sufficiently large so that the gravitational 
perturbation in P2 is negligible. We have a two-body problem between Pq and Pi, where Pq 
is moving with a small velocity a. We can define a weak stability boundary about Pi. This 
can be accomplished by considering a situation where the star Pi is not isolated but included 
within a star cluster. We can think of a new "three-body problem" consisting of Pi, Pq, and 
a third object, B, that mimics the combined gravitational forces from the neighboring stars 
in the cluster. The gravitational perturbation of B together with the gravity of Pi forms a 
weak stability boundary region far from Pi that is located at a finite distance that depends 
on the characteristics of the cluster and the separation of the stars. Pq lies in the transition 
between capture and escape from Pi, when R is sufficiently large and a is sufficiently small. 

5. REMNANT TRANSFER BETWEEN STELLAR SYSTEMS IN OPEN STAR 

CLUSTERS 

Using the framework discussed in the previous sections we now consider the problem 
of remnant transfer between planetary systems in star clusters with low velocity dispersion, 
since low relative velocities are required for the weak capture mechanism. Specifically, we 
consider open clusters which typically have low relative stellar velocities, U ^ 1 km/s. For 
comparison, we note that in older globular clusters, stellar velocity dispersion can reach 
many lO's of km/s. Also for comparison, wc mention that when a cluster starts to disperse, 
the relative distances and relatives velocities between the Sun and neighboring stars increase. 
For example, in the solar neighborhood, the Sun's closest neighbor a-Centauri is 2.6 x 10^ 
AU away (1.28 pc) with a relative velocity of 6 km/s. The latter is significantly higher 
than the ~ 1 km/s required for weak capture, making the transfer of material between the 
two stars very unlikely via weak transfer. Therefore, we will consider relatively young open 
clusters. 

We assume that the stars in the open cluster are approximately uniformly spaced in 
a three-dimensional grid by a distance D and that their distribution is isotropic. Imagine 
a remnant (Pq) in a planetary system {S) passing near the primary planet (P2) having a 
weakly hyperbolic flyby, moving away from the star (Pi) on a pseudo-parabolic trajectory, 
and reaching a distance R — ResdiTii) from the star. At this distance, the small gravitational 
force acting on the remnant due to the central star is roughly comparable with the resultant 
gravitational force from the other stars in the cluster, as described above; as a consequence, 
the motion of the remnant becomes unstable, with small changes to its velocity leading to 
large changes in its trajectory that can either lead to capture or escape from the central 
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stacl. In other words, the sphere of radius R around the central star hes within a weak 
stabihty region and can be thought of as a uniform shce through the more comphcated weak 
stabihty region. 

As shown in Figure [21 due to the structure of the space of parabohc trajectories, for 
any point p on the circle of radius Resc{f^i) around the star Pi, there is a weakly-escaping 
trajectory that will pass by p with velocity a moving away from the star (see also Figure 
1 in Moro-Martm & Malhotra 2005). Assume that the point p lies on the line between the 
two stars. As the remnant (Pq) passes through this point and moves beyond the distance 
Rescij^i) from the star, it goes beyond its weak stability boundary, moving to a region where 
the gravitational force of the star is negligible. The trajectory then continues undisturbed 
until the remnant moves within the weak stability boundary of another star, Pj* in planetary 
system S**, located at a distance Rcapi^l) from Pj*. The remnant can then get captured by 
Pj*, with a periapsis distance that is minimally a collision or maximally a distance Pcap('^i) 
from P^. 

Analogous to planetary system 5", we are assuming that S* has a dominant planet P2* 
at a radial distance A* from the star P^* (see Figure [2]). For weak escape to occur, it is 
necessary that the remnant Pq weakly escapes the star Pi in the same plane of motion 
as that of the dominant planet P2, and that the periapsis distance, Vp, of Pq with respect 
to Pi is approximately the semi-major axis of the dominant planet, A ^ Vp. However, 
when Pq is weakly captured by Pi* and moves to periapsis distance r* with respect to Pi* 
(0 < r* < RcapimD), it need not approach this periapsis within the same plane of motion 
as the dominant planet Pg*. This capture into S* is three-dimensional in nature and the 
remnant Pq can approach the star Pi from any direction. 

The construction of the transfer from the periapsis with respect to Pi (and also P2) to 
the periapsis with respect to P/ can be approximated by determining a solution for Pq in a 
five-body problem between Pq, Pi, P2, Pi and P2*. We can approximate it by piecing together 
the solutions between two three-body problems, Po,Pi,P2 (hereafter TBI), and Po,Pi,P2* 



■^The sensitivity of the motion of Pq at the distance R can be deduced from an analogous four-body 
problem described in Belbruno (2004) and Marsden & Ross (2006) for a transfer to the Moon used by 
the spacecraft Hiten. In this case, the four bodies are the Earth (Pi), the Moon (P2), the Sun (P3) and 
the spacecraft (Pq). The spacecraft leaves the Earth and travels out to roughly 1.5x10^ km where the 
gravitational force of the Sun acting on the spacecraft approximately balances that of the Earth. At this 
location, the motion of the spacecraft is highly sensitive to small differences in velocity, and lies between 
capture and escape from the Earth, i.e. lies at the weak stability boundary between the Earth and Sun. The 
spacecraft then falls back towards the Earth and the Moon entering the weak stability boundary between 
the Moon and the Earth, which finally leads to capture by the Moon. 
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Fig. 2. — Trajectory weakly escaping Pi at a distance R — i?esc("^i); being weakly captured 
by at a distance R — Rcapiml), and moving to periapsis with respect to P* (motion 
projected onto a plane). 

(hereafter TB2). This is done as follows. We know that pseudo-parabolic trajectories in 
the first three-body problem, TBI, exist as described above. In fact, the theory of their 
existence proves there is an infinite number of them, forming the structure of a Cantor set, 
that will travel out to a distance i?ecs(^^i) where the gravitational perturbation due to TB2 
begins to be felt and where the dynamics of Pq becomes sensitive. The velocity of Pq at this 
distance is approximated by the escape velocity from Pi. That is, Pq weakly escapes the 
first system. Let's assume that Pq is then weakly captured by TB2. Its velocity magnitude 
would then approximately be the escape velocity from P^ at the given distance, Rcapi'iTT'i) ■ 
If Pq were moving in the same plane as P2 about Pj*, then by the same theory on pseudo- 
parabolic trajectories, there would exist, by symmetry, an infinite number of trajectories that 
are weakly captured into the system TP2 and fly by P2* in slight hyperbolic state. Then, a 
capture trajectory can be generated in forwards time as an extension of the the weak escape 
trajectory starting at a distance Pcap(^i) from P^*. If, on the other hand, the trajectory 
of Pq does not lie in the same plane of motion of Pg* when it is weakly captured by Pj* , 
then the theoretical results on pseudo-parabolic trajectories cannot be applied. However, 
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the trajectory of Pq can still be extended in forwards time and will move towards P* and fly 
by In this case, it cannot be guaranteed that Pq will fly by Pg* ^ weakly hyperbolic 
manner, although it would be reasonable that this would occur if the time of capture into 
TP2 were sufficiently long. More precisely, one of the following would occur: 

(i) Po remains in a bounded region about P^* without coUision with P*. Then by the Poincare 
recurrance theorem, Pq would eventually fly by Pg* in a shght hyperbolic fashion or coUide 
with P2* 

(ii) Pq collides with P^ 

(iii) Pq is ejected from the P^ — P2*- system. 

The piecing together of solutions at the weak stability boundary regions, for different 
types of trajectories, is discussed in Belbruno 2004 and Marsden & Ross 2006). The piecing 
together occurs at the weak stability boundaries of TBI and T'P2. 

5.1 Location of the Weak Stability Boundary 

We now calculate Pesc(^i) said Rcapiinl) as a function of the stellar mass. For weak 
escape to take place, the velocity, a, of the remnant at the distance Pesc(^i) from the star 
Pi must be sufficiently smaller than the escape velocity at that distance from Pi as well as 
the other stars in the cluster; i.e. it is on the weak stability boundary. At this distance 
the gravitational forces from Pi and from the other stars in the cluster are comparable, and 
the motion of the remnant is unstable and chaotic in nature. Because we are considering 
slow transfer within an open cluster with a characteristic dispersion velocity U ^ 1 km/s, 
we require that a is significantly smaller than [/, i.e. of the order of 0.1 km/s. This is 
much smaller than the nominal values of several km/s used by the Monte Carlo methods in 
previous studies (e.g. Melosh 2003, Adams & Spergel 2005). 

To place the above choice for a in context, we study the velocity distribution of weakly 
escaping test particles from the Solar System, using a three-body problem between the Sun 
(Pi), Jupiter (P2) and a massless particle (Pq). To be consistent with our framework, we 
model this as a planar circular restricted three-body problem, where the test particle moves 
in the same plane of motion as Jupiter, assumed to be in a circular orbit at 5 AU. The 
trajectory of the test particle is numerically integrated by a standard Runge-Kutta scheme 
of order six and numerical accuracy of 10^^ in the scaled coordinates. The initial conditions 
of the test particle is an elliptic trajectory very close to parabolic with periapsis distance 



- 15 - 



rp=5 AU and apoapsis distance ra=40,000 AU. (Note that such orbits are not dissimilar to 
those of known long period comets in the Solar System.) For each numerical integration, 
we assume that Jupiter is at a random point in its orbit when the test particle starts from 
apoapsis at 40,000 AU and falls towards Pi. We record the time at which the test particle 
achieves escape with respect to the Sun (i.e. when the Kepler energy with respect to the 
Sun is positive) after performing a sufficient number of Jupiter fly-bys, and we record the 
resulting hyperbolic excess velocity Voo- Figure [3] shows the distribution of the Voo'. out of 
670 cases, 58% have foo < 0.1 km/s and 79% have v^o < 0.3 km/s. Based on these results 
we will assume the velocity a of the remnant at the distance Resc{mi) from the star to be 
in the range 0.1-0.3 km/s. 




V,„„„,ty (AU) 



Fig. 3. — Velocity distribution of weakly escaping test particles from the Solar System. The 
model is a planar circular restricted three-body problem between the Sun, Jupiter and a 
massless particle. 

For a given a, the location of the weak stability boundary is approximately given by 
Resci^i) = 2Gmi/a'^, where mi is the mass of the star and a is in the range 0.1-0.3 km/s 
(see Figure Hj). Beyond this boundary, we assume that the remnant will move at a constant 
velocity a with respect to the star. 

To allow slow chaotic transfer to a neighboring planetary system, the remnant needs 
to arrive at the distance Rcap{fn\) with a relative velocity with respect to the target star 
{Pi) that is similar or smaller than its parabolic escape velocity at that distance, a* = 
2Gm\/ RcapimX) , where m\ is the mass of the target star; if its velocity is higher than 
a* it will not be captured and and will fly by. Since the relative velocity between stars in 
the cluster is f/ ~ 1 km/s, the remnant that weakly escaped from star Pi moves toward 
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Fig. 4. — Radius of the weak stability boundary as a function of stellar mass. The weak 
escape boundary is defined as Resdjni) = 2Gmi/a'^ with cr = 0.1-0.3 km/s, and the weak 
capture boundary is at Rcapif^i) = 2Gm\/U'^ with U = 1 km/s. The horizontal dashed- 
dotted line indicates the range of cluster sizes, while the dotted line indicates the the mean 
interstellar distance in clusters consisting of N=100 and N=1000 members. 

the target star P* with velocity f/ ± cr (see Figure [2]). Because a is small relative to U 
and can be neglected, the relative velocity of the remnant with respect to the target star is 
PS U. Therefore, weak capture can occur at the distance at which cr* ~ ?7 ~ 1 km/s, i.e. 
Rcap{m\) = 2Gm\/U\ 

Figure m shows Resc{mi) and Rcap{ml) as a function of the stellar mass. The horizontal 
dotted lines indicate the range of cluster sizes and mean interstellar distances for clusters 
consisting of 100 and 1000 members, respectively; such clusters are the birthplaces of a large 
fraction of stars in the Galaxy. The radius of the cluster depends on the number of stars, 
N, and is given by 

/^cluster = lpc(iV/ 100)1/1 (13) 
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-Rciuster is about 2.1 X 10^ and 6.5 x 10^ AU for N=100 and N=1000 members, respectively 
(based on data in Lada & Lada 2003 and Carpenter 2000, and following Adams & Spergel 
2005). We can estimate the average interstellar distance within a cluster, 

D = n"^/^ (14) 

where n = 3N/ {4:TtRI^^^^^^) is the average number density of stars in the cluster. D is about 
7x10^ AU and 10^ AU for a cluster with N=100 and N=1000 members, respectively. 

5.2 Constraints on Stellar Masses for Weak Transfer 

From Figure H] we can set constraints on the stellar mass mi that could allow weak 
escape from Pi to take place. The idea is simple: if for a given a (which as we saw in § 
5.1 is in the range 0.1-0.3 km/s), Resd^i) < D, i.e. the weak stability boundary is located 
within the distance of the next neighboring star P^, then weak transfer is possible because 
at the time the remnant passes near the star Pj*, its velocity is similar to the mean stellar 
velocity dispersion, U (which is a low ~ 1 km/s in open clusters), and there is a significant 
probability of capture (which we quantify in the next section). Conversely, neighbor transfer 
by the process schematically represented in Figure [2] is much less likely to take place if 
Resci^i) > D because at the time the remnant passes near a neighboring star, its velocity 
is too high and as a consequence it will simply fly by. 

Consider a = 0.1 km/s. FigurelHshows that for clusters with 100 members, the condition 
Resc{f^i) < D for weak escape is satisfied for mi < O.4M0, and for clusters with 1000 
members, the condition is satisfied mi < 0.55M©. If we consider the higher but still 
acceptable value a = 0.3 km/s, Figure|l]shows that the stellar mass limits for weak escape are 
mi < 3.5M0 and mi < 5Mq for clusters of 100 members and of 1000 members, respectively. 
A remnant escaping parabolically from any star with larger mass than these limits will 
achieve a velocity of a < 0.1-0.3 km/s only at a distance larger than the mean interstellar 
distance in the cluster and weak transfer is not likely under such conditions. 

Of particular interest is the case of the Sun as the source of the remnants. It has been 
estimated that the Sun's birth cluster consisted of = 2000 ± 1100 members (Adams & 
Laughlin 2001). For a IMq star in such a cluster, we find that the parabolic escape velocity 
at the mean interstellar distance is cr ~ 0.12-0.14 km/s. These values of a certainly lie within 
the range of values of interest for weak escape. We conclude that remnants originating in the 
early Solar System could in principle have met the conditions for weak escape in the Sun's 
birth cluster. 
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5.3 Probability of Weak Capture 

We have established in §5.2 the range of stellar masses that could in principle allow weak 
escape to take place, we now estimate roughly the probability of capture by a neighboring 
planetary system. 

As mentioned previously, because the trajectory of the remnant approaching the target 
system S* is not necessarily in the same plane as that of the orbit of its primary planet 
P2*, and because the capture can be complicated (where Pq could be captured for millions 
of years, moving in a complicated trajectory without planetary collisions), capture is not 
guaranteed even if the renmant Pq falls within the weak stability boundary of the target star 
Pj*. However, because a necessary condition for capture is that the renmant lies within this 
weak stability boundary, we can use the following geometrical considerations to set an upper 
limit for the capture probability. Whether or not the transfer of remnants from one star to 
any star of mass m* takes place will depend on: 

1. The relative capture cross-section of the target star. 

This is given by Cs = Gf{Rcap{ml)/D)'^, where D is the distance between the two stars 
and Reap is described in §5.1 and Figure H] (dashed hne). The factor Gj represents the 
gravitational focusing, given by 

G, = l + (^)^ (15) 

Voo 

where foo is the velocity at infinity and Vesc is the escape velocity at the distance 
Rcap{ml) from P^. 

In the case of weak capture, the term Of increases the cross-section due to enhanced 
gravity focusing. For example, as we saw previously in the case of the Sun and Jupiter, 
f 00 ~ 0.1 — 0.3 km/s, while at the distance Reap = 40, 000 AU, Vesc ~ 0.2 km/s. There- 
fore, in this example, G/ ~ 2, which doubles the capture cross-section. This situation 
occurs in our capture methodology. In the formulation for determining Rcap{ml), the 
remnant has an approximate relative approach velocity to the target star P* of roughly 
U = 1 km/s, which represents the Voo- However, R^ap is determined so that this same 
value of velocity is taken for Vesc from the target star. That is, we are assuming, 
~ Vesc ~ 1 km/s. This implies also that Gf = 2. This is a conservative estimate 
and does not make use of the nature of weak capture dynamics. In this situation, we 
have a trajectory with a Vmf = 1 with respect to Pj*, which goes to a parabolic state 
with respect to P^*. 
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As noted, the value G/ = 2 is conservative. In principle, the value of Gf could be sub- 
stantially increased if at a given value of Rcap{ml), v^o is smaller than the approximate 
value of U = 1 km/s, while the value of Wesc remains the same as U. Dynamically, 
the way to decrease the foo with respect to P^, as the remnant Pq approaches P^, is 
for Pq to decrease its relative velocity. This process has been shown to exist in other 
problems; for example, in the case where we have a spacecraft, Po, transferring from 
the Earth to the Moon on a trajectory that goes to ballistic capture at a given distance 
from the Moon (demonstrated by spacecraft as was discussed previously). Figure 
shows a trajectory to the Moon, going to a periapsis distance of 500 km after 111 
days. When it arrives its velocity is approximately Vesc- However, its v^o goes from a 
value of 1 km/s to 0. This is reflected in Figure [6] of the Kepler energy Km of Pq with 
respect to the Moon along this transfer. At a sufficiently far distance from the Moon, 
where Voo ~ V^Ke, Voo approaches zero. In the future, we would like to examine the 
distribution of the f oo in the WSB about P^ and thereby get a better understanding 
of the Gf and the enhanced capture probabilities. 

# 
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Fig. 5. — Trajectory from the Earth to the Moon going to the lunar WSB. 

2. The number of potential targets. This depends on the probability of finding a star 
of a given mass in the cluster; we call this PiMpij^i) which describes the initial 
mass function (IMF) of the cluster. Observations of many different star clusters (large 
clusters like the Trapezium, smaller cluster like Taurus and even older field stars) find 
very similar IMFs (down to the Hydrogen burning limit at ~ 0.1 Mq; Lada & Lada 
2003 and refereces therein). To calculate PiMpif^i) we adopt the IMF of the Trapezium 
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Fig. 6. — Caption: Variation of the Kepler energy, Km- 

cluster from Lada & Lada (2003), which is characterized by a broken power-law, given 
by ^(M) = 6M-2•^ from 0.6-20 M© and ^(M) = ^M'^'S from 0.1-0.6 Mq, where 
(^(M)dM is the number of stars with mass (M, M+dM). [There is a steep decline into 
the substellar brown dwarf regime and a possible second peak but we will ignore objects 
below the hydrogen burning limit]. To calculate Pjmf{^i) (square symbols in Figure 
[7]), we use a logarithmic binning of masses (with d(logM)=0.1), and normalize the 
distribution to unity, which gives ^i=0.19 and ^2=0.34. 

An upper limit to the probability that a remnant escaping from a star of mass mi (within 
the range described in §5.1) will get captured by a neighboring star of mass m* is approxi- 
mately given by 2(_Rcap("^i)/-Dm*-mi)^) where -Dm* -mi is the average distance between two 
stars of masses ml and mi, respectively. The simplest case is when both stars have equal 
masses mi = mi. In this case, the average interestellar distance would be ~ {^/n-miY^^, 
where n^i the average number density of stars with mass mi, rimi = A^mi/(4/3)7ri?^;„g^g^, and 
is the total number of stars in the cluster with mass mi, A^^i = N x PiMpi'n^i)- Because 
= nx PjMFirrii), "we get that D„ii ~ -D x (P/a/f("^i)) This means that the transfer 
probability between two stars of equal mass is given by 2{Rcap{rrL[) / DY x (P/Mi?(mJ))^/^, 
where D is the average distance between any two stars in the cluster (regardless of their 
mass). The resulting capture probability is shown in Figure [7] as a solid line, with values 
ranging from 10~^ to 10~^. The capture probability between two planetary systems with 
solar-type central stars (MJ = Mi = 1 Mq) is 8.1x10"^ and 1.7x10"^ for a cluster of 1000 
and 100 members, respectively. 
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Fig. 7. — For a given source of remnants, the probability of capture by a neighbor 

of equal mass, shown in solid line, is given by 2{Rcap{fn*i) / DY x (P/mf(^^i))^^^; where 
2{Rcap{m\) / DY is shown in dotted hne and PiMF{jn\) is the Trapezium IMF normalized to 
unity, shown as square symbols (Lada & Lada 2003). 



5.4 Estimate of the Number of Weak Transfer Events 



To calculate the total number of remnants that could get transferred between two neigh- 
boring planetary systems, the capture probability in §5.3 needs to be multiplied by the 
number of remnants, Nr, that a planetary system may eject before the cluster disperses. 
The main uncertainties in assessing whether weak transfer is a viable method for remnant 
transfer lie in estimating this number. Given the large uncertainties, in this sub-section (and 
the remainder of this paper) we will only estimate Nr for our Solar System, as it is the 
only planetary system for which observations and dynamical models enable us to make an 
educated estimate. 

An estimate of was given in Adams & Spergel (2005): assuming that the early solar 
nebula contained about 50-100 of heavy elements, from which about one third was left- 
over from the planet formation process, and about one third of this left-over material was 
ejected from the planetary system during encounters with the giant planets, they estimated 
that during the time of planet formation, which lasted ~ 10 Myr, at least Mji ~ 1 of 
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rocky material could have been ejected from a planetary system. Because clusters remain 
bound during 10-100 Myr, some of this rocky material could in principle get captured by 
other stars in the cluster. Of particular interst are the remnants ~10 kg, large enough to 
shield any potential biological material from the hazards of radiation in deep space and 
from the impact on the surface of a terrestrial planet (Horneck 1993; Nicholson et al. 2000; 
Benardini et al. 2003; Melosh 2003). 

5.4.1. Estimating from the Oort Cloud 

An estimate of Nji can be derived from observations and dynamical models of the Oort 
Cloud comets, which are weakly bound to the Solar System and are therefore representative 
of a population of remnants that may have been located on its weak stability boundary, 
subject to weak escape. 

Oort Cloud formation scenario: 

In a recent paper, Brasser, Duncan & Levison (2006) proposed the following scenario for 
Oort Cloud formation. Before the gas in the solar protoplanetary nebula dispersed, the gas 
giant planets Jupiter and Saturn formed, and subsequently scattered the Jupiter-Saturn zone 
planetesimals out to large distances. This process happened relatively quickly, before the 
Sun left its maternal stellar cluster and the cluster gas dispersed. At these distances, the 
planetesimals were subject to the gravitational perturbations of the cluster gas and stars; 
these perturbations cause a Kozai-hke effect of coupled eccentricity-inchnation oscillations. 
These long period oscillations in a slowly changing gravitational potential of the cluster, 
resulted in the lifting of the pericenter of the planetesimals' orbits beyond the orbit of 
Saturn. Close encounters with stars, in particular near the center of the cluster, also had 
the effect of increasing the pcriccntcrs. Planetesimals that achieved pcriccnters ^ 10 AU 
were safe from complete ejection from the Solar System, and subsequently formed the Oort 
Cloud. (This model does not consider the ice giants Uranus and Neptune because they 
hkely formed after the solar protoplanetary gas dispersed (>10 Myr), and probably after the 
Sun left the cluster.) Brasser et al. best estimate cluster model that correctly explains the 
current orbits of Sedna and 2000 CR105 and that makes Sedna a typical Oort Cloud member 
consists on a stellar cluster of 284 members, with a total mass in stars and gas of 419 Mq, a 
velocity dispersion of 1.755 km/s, a crossing time of 0.056 Myr, a Plummer radius of 0.1 pc 
and a tidal radius of the Sun of 0.013pc. In this model, 2-18% of the planetesimals in the 
Jupiter-Saturn region became part of the primordial Oort Cloud. 

We adopt the above scenario to calculate Nr. First, we estimate the number of plan- 
etesimals in the Jupiter-Saturn zone, which we assume was the 4-12 AU heliocentric distance 
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zone; we assume that 2-18% of those ended up in the Oort Cloud. 

Total mass of solids in the primordial 4-12 AU region: 

To calculate how many planetesimals formed in this region, we estimate the total mass 
in solids (i.e. excluding H and He) that could have been available for the formation of 
planetesimals in the 4-12 AU region. We adopt the minimum mass solar nebula (MMSN), 
the minimum mass that the solar protoplanetary disk must have contained in order to form 
all the planets in the Solar System, given by the surface density E = Eo(a/40^C/)~^/^, 
where the surface density of dust (i.e. sohds) has the coefficient Eq = Eq^ = 0.1 g/cm^, 
(Weidenschilling 1977, Hayashi 1981). Integrating between 4 and 12 AU, we find the total 
mass in sohds, 10^^ g. 

Planetesimal size distribution: 

With the above estimate for the total mass in solids, we can calculate the number of plan- 
etesimals by adopting a planetesimal size distribution function representative of the early 
Solar System. This size distribution is highly uncertain, but can be constrained roughly 
from observations and coagulation models. 

The current best observational estimate for the size distribution of outer Solar System 
planetesimals is for the trans-Neptunian Kuiper belt bodies studied by Bernstein et al. (2004) . 
Theoretical estimates are based on planetesimal coagulation models, reviewed in Kenyon et 
al. (2007). Here we summarize these results, and use them to guide our estimates for Nr. 

• Observations of Kuiper Belt bodies show, broadly, two dynamical classes, the Classical 
Kuiper Belt (CKB) with low inclination, low eccentricity orbits and the Excited Kuiper 
Belt (EKB) with modcratc-to-high orbital inclinations and eccentricities. Bernstein et 
al. (2004) find that the size distribution functions of these two classes arc different at a 
96% confidence level. Each class is well fitted by a 'broken' power law, dN/ dD oc D~'^i 
for D > Dq and dN/dD oc D"^^ if D<Do, of different power law indices at the small 
and large sizes: ^^>5.85 and q^^^^2.9 for the CKB, and gf ^^Ri4.3 and gf ^^<2.8 
for the EKB. The break in the power laws is at £)o ~ 100 km. 

• The CKB has fewer large objects than the EKB; the maximum planetesimal size in the 
CKB is 60 times smaller than the maximum planetesimal size in the EKB. The EKB 
contains fewer small objects than the CKB. This differences indicate that the Excited 
population is likely "older - has undergone more coUisional accretion and erosion - 
than the Classical population. This suggests that the EKB objects may have formed 
at smaller heliocentric distances than those in the CKB. 
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• For both populations, the size distribution of the small bodies is shallower than ex- 
pected in a coUisional cascade, q2 < 3.5, where the latter value is a theoretical estimate 
for collisional cascades (Dohnanyi 1969). The break to a shallower size distribution 
occurs at D < 100 km, a size range susceptible to collisional destruction (Pan & Sari 
2005). The observed strong depletion of the small bodies, and the fact that the col- 
lisional lifetimes in the present-day Kuiper belt are longer than the age of the Solar 
System, indicates that in the past both the CKB and EKB were more massive and 
richer in small bodies and that the present state is the result of an advanced erosional 
process. 

• Theoretical models of planetesimal coagulation in the outer Solar System also find a 
broken power law size distribution function, with parameter values qi ~ 2.7-3.3 and 
q2 — 3.5, and a break diameter Do ~ 1km (reviewed in Kenyon et al. 2007). These do 
not match the observations of the present-day Kuiper belt, but may be representative 
of its very early stage during the planet formation era. 

Based on the above studies, we adopt three size distribution functions for our calculation 

• Case A: qi = 4.3, q2 = 3.5, Dq — 100 km, Dj^ax — 2000 km (~ Pluto's size), D^j„ = 
1 fim (~ dust blow-out size). 

This distribution has the power law index of a collisional cascade at the small size end, 
and that of the EKB at the large size end, with a break diameter consistent with that 
of the present-day Kuiper belt. The rough scenario this case reflects is as follows: In 
the Jupiter-Saturn zone, the accretion of large planetesimals proceeded to make the 
large-size end similar to that found in the present-day EKB, whereas at the small size 
end, the dynamical stirring by the large bodies produced a classical collisional cascade. 

• Cases B and C qi = 3.3 and 2.7, q2 = 3.5, Dq = 2 km, T>rnax = 2000 km (~ Pluto's 
size), I^rnin = 1 Aini (~ dust blow-out size): 

These size distributions are derived from theoretical coagulation models. 

• Case D: qi = 4.3, q2 = 1.1, -Do = 100 km, D^ax = 2000 km (~ Pluto's size), Dmin — 
1 fim (~ dust blow-out size): 

This size distribution represents perhaps the worst case scenario in which the depletion 
of the small bodies took place very early on, before the objects were thrown into the 
Oort Cloud. The parameters here are within the range that Bernstein et al. (2004) 
find for the EKB. This size distribution is also similar to that found in models for the 
primordial asteroid belt (gi = 4.5, q2 — 1.2, Dq ~ 100 km, Bottke et al. 2005). 
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Number of remnants with masses > 10 kg that could have been subject to weak escape: 
With the above estimates of the total mass in sohds and of the size distributions, we can 
finally calculate the number of planetesimals with masses > 10 kg, equivalently, diameter 
D > 26 cm (assuming p = 1 g/cm^^), that would have existed in the 4-12 AU region 
at the time of Jupiter and Saturn formation. The results are shown in Tabic 1. The 
scenario described in Brasser et al. (2006) finds that 2-18% of these planetesimals would 
have tranferred to the primordial Oort Cloud; these are also the bodies potentially available 
for weak escape from the early Solar System, hence they provide our estimate for Nr, also 
given in Table 1 for each of our cases. We find that Nr is in the range 10^^-10^°'^ for Cases A, 
B and C (which are based on a coUisional cascade for the small objects), but is only 10^-10^ 
for Case D (which has a much shallower size distribution for the small bodies, similar to the 
observed Kuiper belt). It is clear that Nr is very sensitive to the size distribution function, 
i.e. power law index used. 
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5.4.2 Number of Weak Transfer Events in the Early Solar System 



Now we estimate the number of weak transfer events from the early Solar System to the 
nearest solar-type star in the cluster (assuming it also harbored a planetary system) by 
multiplying Nn with the capture probabihty calculated in section 5.3 (see Fig. 5). For = 
Ml = 1 Mq, the capture probability is 8.1 xlO^^ for N=1000 and l.TxlQ-^ for N=100. The 
result is also listed in Table 1. The number of weak transfer events between two solar-type 
stars for Cases A, B and C are in the range lO^^-^-lO^^-^, and for Case D is 900-17000. As 
mentioned before, because the location of the renmant inside the weak stability boundary 
does not guarantee capture, these estimates should be regarded as upper limits. 



5.5 Implications for Lithopanspermia 



Because planetesimals could potentially harbor the chemical compounds that constitute 
the building blocks of hfe, it is of interest that the results above indicate that these mate- 
rials could have been transferred in significant or even large quantities between the Solar 
System and other solar-type stars in its maternal cluster. For the case of binary systems, 
Adams & Spergel (2005) found that clusters of N = 30-1000 could experience billions to 
trillions of capture events among their binary members. However, to consider the transfer 
of microorganisms that developed in the Solar System, a necessary condition would be that 
life could develop during the 10-100 Myr that the cluster remained bounded. Radiomet- 
ric measurements of hafnium and tungsten isotopes in meteorites indicate that the bulk of 
the metal-sihcate separation in the Solar System occurred within the first 30 Myr of the 
solar nebula hfetime, with most of the Earth's core accreted during the first 10 Myr (Yin 
et al. 2002; Kleine et al. 2002). It is not known when the conditions for hfe to develop on 
Earth were met. There is evidence, albeit sparse, that basic habitable conditions (existence 
of continents and prevalance of liquid water) were established on Earth within ~ 150 Myr of 
Solar System formation (e.g., Catcs & Mojzsis 2007). The oldest Earth rocks are ~ 3.9 Gyr 
old and some authors claim that the first traces of life appear just 100 Myr later (Mojzsis 
et al. 1996), but this result is controversial (Fedo & Whitehouse 2002). Similarly, there is 
controversy surrounding the finding of modern cyanobacteria in rocks 3.465 Gyr old (Schopf 
1993; Brasier et al. 2002); undisputed is evidence of cyanobacteria and primitive eucaryotes 
in rocks 2.7 Gyr old (Brocks et al. 1999). 



Transfer Timescales 
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For the transfer of microorganisms from the Solar System to other systems, we need to 
consider not only the timescale for life to develop, but also the timescale for weak transfer 
and how it compares to the timescale for microorganism survival in deep space. 

• Timescale for Ejection: The time, Tr, for a remnant to exit the Solar System, i.e. to 
move from its periapsis (r^) with respect to the central star to the distance Resdjni) 
can be is estimated using Barker's equation which yields the time of flight along a 
parabolic trajectory, 

where 

v p 
L = ^tan( — ), u = arccos(— 1). 

2 Resc 

The variable u is the true anomaly, and p = 2rp{= 2A) is the semi-latus rectum. For 
p <^ Resc and using trigonometric identi ties it can be shown that L varies approximately 
independently from p : L = ^2 -Resc a /l — ^r^— ~ V^Resc ■ Therefore, the variation 

y ZiXesc 

of Tr as a function of is negligible, so without loss of generality, we can examine 
the variation of Tr as a function of Resc for a fixed periapsis distance of Vp ^ 5 AU: 
for stellar masses of 0.1 Mq, 1 Mq and 10 Mq, Figure H] shows that Resc = 1.8x10^ 
AU, 1.8x10^ AU and 1.8x10^ AU, respectively, which yields Tr ^ 0.6, 6 and 60 Myr, 
respectively. 

• Timescale for interstellar transfer. For a remnant moving at the low velocity of 0.1 
km/s, as required by the weak transfer mechanism, it will take about 3.3-4.8 Myr to 
reach a neighboring star located at D ^ 7x10^-10^ AU. 

• Timescale to land on a terrestrial planet: Because it is not required that the captured 
remnant approaches the star in the ecliptic plane, multiple periapsis passages about 
the target star would be needed before the remnant can collide with a planet. This 
can take of the order of tens of millions of years. 

• Timescale for erosion from interplanetary dust: In the interplanetary dust environment 
found in the Solar System, a meter-size rock will be ablated and eroded on a timescale 
of about 0.02-0.23 Myr, short compared to the time it takes to be ejected from the 
system (about 4 Myr; Napier 2004). It is likely that the early Solar System was 
significantly more dusty than it is today, thus survival times of remnants would have 
been even shorter, and much less than the several Myr transfer timescales. 



Even if we were to assume that life arose on Earth before the stellar cluster dispersed, 
given the ejection and transfer timescales described above and the survival timescale of 
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dormant microoganisms, it seems unlikely that microorganisms could have been transferred 
even among the closest neighbors in the cluster via the weak transfer mechanism described 
here. However, our results that significant quantities of solid material could have been 
transferred via weak transfer between nearest neighbors in the Sun's birth cluster, are worth 
further investigation for the exchange of chemical compounds that constitute the building 
blocks of life and which could potentially have significantly longer survival lifetimes. 

6. CONCLUSIONS 

Could life on Earth have been transferred to other planetary systems within the first 
few Myr of the Solar System evolution, when the Sun was still embedded in its maternal 
aggregate?; and vice versa, could life on Earth have been originated beyond the boundaries 
of our Solar System? 

In this paper we have described a dynamical mechanism that yields very low velocity 
chaotic escape of remnants from a planetary system using parabolic trajectories, and its 
reverse process of chaotic capture. These two processes provide a mechanism for minimal 
energy transfer of remnants between planetary systems. We have applied this mechanism 
to the problem of planetesimal transfer between planetary systems in an open star cluster, 
where the relative velocities between the stars are sufficiently low (~ 1 km/s) to allow slow 
escape and capture. Based on geometrical considerations, we have estimated upper limits 
to the probability of transfer that depend on the cluster properties and the masses of the 
source and target stars. Using these probabilities, adopting the Oort Cloud formation models 
from Brasser et al. (2006), and adopting a range of planetesimal size distributions derived 
from observations and theoretical models (Bernstein et al. 2004 and Kenyon et al. 2007, 
respectively), we estimate the number of weak transfer events from the early Solar System 
to the nearest solar-type star in the cluster (assuming it also harbored a planetary system). 
This estimate is most sensitive to the power law index of the size distribution of small 
bodies: the number of weak transfer events could be could be as large as lO^^-^-lO^^-^ if the 
size distribution of the small bodies follow a classical coUisional cascade, or as small as 900- 
17000 if we adopt a shallower power law index for the size distribution of the small bodies. 
To determine whether the weak transfer process described in this paper could have been a 
viable mechanism for the transfer of remnants between the Solar System and other stars in 
the cluster, further progress needs to be made in the understanding of the dynamical and 
coUisional history of the early Solar System. 

Comets and asteroids have been suggested to be a possible source of the chemical com- 
pounds that constitute the building blocks of life on Earth. The results summarized above 
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are therefore of interest because they indicate that there is the possibihty (depending on the 
shape of the size distribution) that planetesimal material could have been transferred via 
the weak transfer mechanism in large quantities between the Solar System and other stars 
in its maternal cluster. 
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Table 1. Estimated Number of Weak Transfer Events 



^For all the 

^^D>26cm is 



Planetesimal Size ^D>26cm^ ^wte"^ ^wte"^ 



=^Wkg if p 



Distribution'^ (N=100) (N=1000) ^^-^ ' 
primo rdial Oo 

A: gi = 4.3, 02 = 3.5 1.8x10^^ 3.6x10^^^-3.2x102° 6.1x10^^-5.5x10^^ 2.9 x 10^^-2.6 xlOi^^, 

r^WTE IS th 

Do = 100 km J i\T 1 nnn 

and N=1000 ] 

B: = 3.3, 92 = 3.5 2.7x10^0 5.4x10^^.9x10^9 9.2x10^^-8.3x10^^ 4.4 x 10^^-3.9 x J(^^ ^q_4 ^ 
Do = 2 km 

C: gi = 2.7, g2 = 3.5 8.1x10^^ 1.6x10^^-1.4x10^^ 2.7x10^^-2.5x10^^ 1.3xl0^3_i_2xl0^^ 
Do = 2 km 

D: gi = 4.3, 92 = 1.1 5.5x10^ 1.1x10^-9.9x10^ 1.9x10^-1.7x10^ 8.9x10^-8.0x10=^ 
Do = 100 km 



